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Abstract 

We compute the Bohmian trajectories of the incoming scattering plane waves for 
Klein's potential step in explicit form. For finite norm incoming scattering solutions we 
derive their asymptotic space-time localization and we compute some Bohmian trajecto- 
ries numerically. The paradox, which appears in the traditional treatments of the problem 
based on the outgoing scattering asymptotics, is absent. 

PACS: 03.65.Bz, 03.65.Pm 

1 Introduction 



Bohmian mechanics ]lj []2| attempts to reconcile quantum mechanics first with the notion of 
observation independent properties of physical systems ("realism") and second with strictly 
deterministic basic laws connecting these properties ("determinism"). While in standard quan- 
tum mechanics a system acquires a specific property only upon measurement by random quan- 
tum jumps, within Bohmian mechanics each individual system has all its possible properties 
completely specified independently of any measurement. This is achieved by completing the 
quantum states through "hidden variables" , which determine all outcomes of experiments with 
individual systems. Whether the laws, which connect hidden variables, can be tested exper- 
imentally, is still under debate. In any case the basic aim of Bohmian mechanics has 
been appreciated by a growing community, [|| and apart from all fundamental controversies, 
Bohmian trajectories undoubtedly help to visualize the time evolution of wave functions. 

While the main bulk of work within Bohmian mechanics has been devoted to the Schrodinger 
dynamics, there exist some results concerning relativistic quantum mechanics. Among these is 
Holland's work on the Klein paradox. Q In this work the Bohmian trajectories are computed, 
which follow from the treatment of Klein's Paradox as it is given by Bjorken and Drell. 
Bjorken and Drell present a plane wave solution of the Dirac equation for an electron, which 
is exposed to a sufficiently high, one dimensional potential step. They talk about the vari- 
ous parts of their solution as incoming, transmitted, and reflected waves. The paradox then 
arises because the reflected probability current is greater than the incoming current and the 
transition current is directed towards the potential step. Holland Q determined the Bohmian 
trajectories associated with the plane wave solution of Bjorken and Drell. He showed that they 
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are time like, future directed and do not intersect. Furthermore they do not begin or end at 
finite time. Thus there is no indication for the production or annihilation of electron-positron 
pairs in these trajectories. However, there is something strange with them. They emerge from 
the region with nonzero potential instead of moving into this region in the course of time. 
Furthermore the trajectories of the physically interpretable finite norm wave packets have not 
been considered at all in ||. 

Already long ago Bongaarts and Ruijsenaars || pointed out the reason for the paradoxically 
directed transition current of Bjorken and Drell's treatment: the chosen sign for the momentum 
parameter in the region with nonzero potential. As a consequence, the plane waves of Bjorken 
and Drell obey the boundary condition of an outgoing scattering solution instead of an incoming 
one: A wave packet superposition of these waves, is localized within the potential free half space 
and propagates away from the potential step at large positive times. At large negative times 
the wave packet approaches the potential step from both sides. Thus the surrounding talk of 
"incoming", "transmitted", and "reflected" waves, which was used by Bjorken and Drell, is 
mathematically unjustified and misleading. 

In contrast to these outgoing wave packets, a wave packet superposition of the corresponding 
incoming plane wave solutions, is localized within the potential free half space and propagates 
towards the potential step at large negative times. At large positive times it moves away from 
the potential step in both directions. Thus in this case the picture of an incoming, transmitted, 
and reflected wave is based on mathematical facts. Elaborate numerical studies of the evolution 
of Gaussian wave packets have been made in order to confirm the completely different behaviour 
of the two types of wave packets. Q As a result, the paradox of a negative transition current, 
as formulated by Bjorken and Drell and investigated by Holland, disappears, if the incoming 
scattering solutions are used in place of the outgoing ones. 

In this work we first discuss the Bohmian mechanics, associated with the incoming scattering 
solutions of Klein's paradox, which initially move in the potential free region towards the 
potential step. For the building block plane waves we prove that all Bohmian trajectories move 
into the region with nonzero potential. Second, we consider the finite norm wave packets. 
We derive their asymptotic space-time localization properties and compute some trajectories 
numerically. We find that at large negative times, all trajectories, which carry a substantial 
part of the total norm, are located in the region with zero potential and are directed towards 
the potential step, where some of them are reflected and some are transmitted. At large positive 
times these trajectories move away from the potential step. No indication for pair creation at 
the step can be found. There remains, however, as will be shown here, a physically questionable 
acceleration of slow packets upon passing the potential step. 

2 The 2d Dirac equation 

The set of space time points is assumed to be M := K 2 . Let id,M =■ (a: , a; 1 ) denote the 
standard chart of M. The associated tangent frame is d := (do,di). A Minkowskian metric on 
M is defined by G p (d M , c\,) :— 5®5® — 5^51 for each p 6 M. Let e = (ei, e2) be an orthonormal 
basis of a 2dimensional (2-d) C- vector space W with scalar product S : W x W — > C. The linear 
mappings 7 M : W — > W are defined by 7°(eo) = ei, 7°(ei) = e and 7 1 (e ) = ei, 7 1 (ei) = — eo- 
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A diffcrentiable function ip : M — ► W with 

(d„ - i^A^ tf, = Kip (onM) (1) 

is called a classical solution of the Dirac equation with continuous external potential A = 
A^dx^ = {Afj) and Compton length kT 1 £ R>o- The differential geometric notation is as in 

Let ip 1 and ip 2 : M — * C denote the (spinor-) component functions of the function ip : 
M — > W with respect to the basis e, i.e. ip —: eiip 1 + t2ip 2 ■ The (Lorentz invariant) indefinite 
inner product L of the spinor space W in terms of its scalar product S reads L(v,w) := 
S [v, 7°w) = v 2 *w x + v x *w 2 . The current vector field j^p := j°<9o + j^di = of a function 

ip : M W is defined by ^ := L(^,7^)^ = (l^ 1 ) 2 + I^T) <9 + (|V^| 2 - l^ 2 | 2 ) 9 i- Duc 

to G(jip,ji/,) — 4|-0 1 t/' 2 | 2 > 0, the current is nowhere space-like. Where is nonzero, it is 
future oriented. If tp is a classical solution of the Dirac equation (0), then div(j^) = 0. 

For sufficiently regular potential, Q from div(j^) = and j° > a probability structure 
for the spaces E T := {p € M : a; (p) = r} with r€l can be established as follows. The norm 
of the restriction tp T of a classical solution ip to S T for any r € K is defined by 



For solutions with ||^o|| < °o the equation HVvll = HV'oll holds for any r G R. Therefore the 
density 



W.r := I, , ..a with £ : = 2:1 

11-00 1| 



is a probability density on E T . 

It has been suggested in sect. 12.2 of Q, || that, due to div(j^) = 0, the density p^ }T 
is the transport of p^,o from So to E r along the flow lines of j^. The phenomenon may be 
visualised as the evolution of the mass distribution of a cloud of dust along the individual 
particle trajectories. In consequence, the set of flow lines of j^,, the Bohmian trajectories, have 
been taken seriously as the possible particle world lines, i.e. each (one particle) system in the 
quantum state represented by ip is supposed to realise one of the flow lines of in the course 
of time. 

A general set of potentials and initial conditions (including singular ones) seems to be 
unknown such that the global Bohmian trajectories densely fibre supp(j^,). (A global trajectory 
is one, which extends both unboundedly into the past and into the future. A dense fibration 
is such that the set X C supp(j^) of points, which lie on a global Bohmian trajectory, obeys 
/°V>.r(A n £ T ) = 1 for any r € R.) The analogous problem in the Schrodinger case has been 
described by Berndl in || and resolved in |Q for a wide class of potentials. 

Lacking such general results for the Dirac equation we confine ourselves to the very specific 
case of a discontinuous potential step. Let 0<^el and 6 : M — > R denote the step function 
Q(x > 0) := 1 and Q(x < 0) := 0. Then the differential 1-form A := A dx° with -jfeAo := 
V ■ (O o x 1 ), defined on U := \p : x 1 {p) ^ 0}, is introduced as an external potential into the 
restriction of (m) to the domain U. Only classical solutions of the restricted equation with a 
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continuous extension to M are taken into consideration. This yields the following system of 
partial differential equations for differentiable component functions ip l :[/—*■ C with continuous 
extension to M. 

ld ° { v> 2 ) = { K id 1 + ve&) ) U 2 ) (on u) (2) 

For any such solution ip, the current jj, is continuous on M and differentiable on U. On U 
the continuity equation div(j^) = holds. 

3 Localization of free wave packets 

The finite norm solutions of (||) with V > will be constructed from the plane wave solutions 
of (Fy) in the case A = 0. Thus this case is summarised first. See e.g. ]f3[. 



N otation 1 Let uJ : R -> R >0 urai/i uJ(fc) := Vk 2 + fc 2 and fi : R -> R >0 w«£/i fi(fe) := 
A/w(fc) + fc. T/ien it:R-vW>:R-*-W with 

Remark 2 For ^ ooi/i (u(k),u(—k)) and (v(k),v{—k)) is a basis of W . The following 
relations hold. 

S(u(k), u(k)) = 2ZJ(k), S(v(k),v(k)) = 257(fc), 

S(u(k),u(—k)) = 2k, S(w(fc), «(-*;)) = 2k, 

L(u(fc), 7 1 u(fc)) = 2fc, £(«(*), 7 X «(A)) = 2*, 

L^fe)^ 1 ^-*;)) = 0, LKfej^V-Zc)) = 0. 

Lemma 3 Let k £ R \ and Ze£ </ie function f : M ^ W be differentiable. Then a) and b) 
hold. 

a) exjp(—iuj(k)x°)f(x 1 ) : M —* W solves the Dirac equation ^) with A = if and only if 
for some a, [3 G C 

/ (a: 1 ) = a exp(ikx )u(k) + /3exp(— ikx 1 )^— k). 

b) exp(iUJ(k)x°) f (a; 1 ) : M — > W solves the free Dirac equation ^) with A — if and only 
if for some a, f3 € C 

/ (a; 1 ) = a exp(— ikx 1 )v(k) + (3 exp(ikx 1 )v(~k). 



Notation 4 For k S R we denote \ph^Uu := exp(— i(uj(k)x° — fcx 1 )u(fc), and y/2wVk := 
exp(i(ZU(fc)x° — kx 1 )v(k). 

The frequency of these plane wave solutions belongs to (— oo, — k)U(k, oo). To each frequency 
within this range a 2-dimensional subspace of single frequency solutions to (|f|) exists. The space 
of solutions with frequency ZU(fc) > k is spanned by ([/&, U-k) and the space of solutions with 
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frequency —ZJ(k) < — k is spanned by (Vk, V-k)- Both Uk and Vk are constant along the space- 
like phase velocity vector field / := <9 + ^jp-di. The current vector field, associated with both 
Uk and T4, is future oriented, time like, and constant: 

3u h = jv k = -(uJ(k)d + kdi). 

TT 

From the plane wave solutions Uk and Vk of (|l|) finite norm wave packets are formed by 
superposition. Let / C R be a compact interval and let the function a : / — > C be continuous. 
Then the functions from I x M into W, with either (k,p) i— > a(k)Uk(p) or (fc,p) i— » a(A:)T4 (p), 
first are continuous and second have continuous partial derivatives with respect to x^. There- 
fore the functions from M into W with either p L a(k)Uk(p)dk orpn fja(k)Vk(p)dk 
are differentiable and the differentiation may be interchanged with the integration and these 
functions are solutions of ([!]) with A = 0. 

Notation 5 For a continuous function a : I — > C, defined on a compact real interval I, wave 
packet solutions U [a] and V [a] of ^) are defined by 

ft dk 
U [a] := / d/i(k)a(k)Uk and Via] := / d[i(k)a(k)Vk with d/i(k) :— — 
Ji Ji 2w(fc) 

U [a] is called positive frequency packet, V [a] is called negative frequency packet. 

The movement of "narrow" wave packets U [a] and V [a] can be made plausible by replacing 
the function ZJ in the factor exp(— ix°oJ) by its tangent approximation at a point fco from the 
domain /. This yields U [a] ~ P+ [a] U ko and V [a] ~ P_ [a] V ko , with w° := and 

P± [a] := dn(k)a(k) exp [qp*(Ai - fco) • (^a; - a; 1 )] • 

The functions P± [a], which modulate the plane waves Uk and I4 , are constant along the 
(future directed, time like) group velocity vector field g := <9o + w° • d\ on M. Thus the sign of 
fco determines the direction of propagation of P± [a] and in case of fco > the approximations 
to both U [a] and V [a] propagate towards growing x 1 . 

A more conclusive derivation of the space-time localization of the wave packets U [a] and 
V [a] follows from proposition (3.1) of ref. 0. In the present case of one space dimension this 
proposition reads as follows. 

Proposition 6 Let T : L 2 (R) — > L 2 (R) denote the Fourier transformation, formally given by 

1 f°° 

J- [f) : fc i— > / exp(— ikx)f(x)dx. 



/2tt 

Let uj : R — > R be twice continuously differentiable. The first derivative of lu is denoted as 
U! 1 . Define for any i £ R the unitary time evolution operator Ut '■ L 2 (M) — > L 2 (R) through 
ut{f) '■= ^ r_1 (exp(— iujt)T(f)). Let f S L 2 (M) and v\,V2 £ R be such that v\ < u/(fc) < vi for 
all fc € supp.F(/). Then 

/tv\ /"OO 
\u t (f)(x)\ 2 dx = lim / | Ut (/)(x)| 2 da; = 0. 
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Remark 7 Since ||/||^2 = ||ut(/)||x,3 f or oil ' £ R, the above statement is equivalent to 



lim 



\u t (f)(x)\ 2 dx = \\ff 



L 2 



u t(f) localises for t — > oo within the interval t ■ [v\, vq\. 

Remark 8 The localization of Ut(f) for t — > — oo can be obtained from the limit t — » oo of the 
evolution u t := U-t> which has the frequency function uj :— In this case — V2 < uj'(k) < —V\. 
Thus 



lim 

t — KX) 



— tvi 



\ u t(f)(x)\ dx 



IV2 



lim 

t — > — oo 



\ u t{f){x)\ dx 



tv 2 



follows. Both limits are covered by 

WfWl- 



lim 

i— »±oc 



\ u t{f)(x)\ dx. 



>t-[v 1 ,v 2 ] 

An application of proposition (Q) to the component functions of U [a] and V [a] yields the 
following localization for p^^ T . 

Proposition 9 v g := to' : M — > M, k i— > =^y denotes the group velocity of 57. The function 
a : [kifk^] — > C be continuous. The real numbers v\ and v^ are chosen such that vi < v g (ki) 
and v g {k2) < f2- For ip G {U [a] , V [a]} then 



lim 

T — >±QO 



t[vi,v 2 ] 



1. 



Proof. We shall check first that the assumptions of proposition (||) hold for the component 
functions of our wave packets. Consider the positive frequency case. Since the functions and 
uj are continuous and since = is bounded, the function 



q(fc)n(-(-l) i fc) 
2U(k) 





for k\ < k < &2 
otherwise 



belongs to L 2 (R). Now U [a] 1 is given by U [a}\ — T 1 (exp(— iujr)g l ) . As uj has continuous 
derivatives of arbitrary order, the evolution is of the type of proposition (^|). The second 
derivative W"{k) = =fep > implies that the group velocity function k > uj'(k) — is 
strictly increasing. From this one obtains the bounds v\ < U'(k) < V2 for any k G [fci, A^]. Thus 
proposition (||) yields 

2 



U\a\ 



I St 



L 2 
2 



lim 

T — >±00 



T-[Vl,V 2 ] 

2 



\U[a]\T,x)\ dx. 



From this and ||{7[a] r || = ||[/[a]*|| L2 + ||C/[a]£ || 2 the statement follows. The case of negative 
frequency packets is analogous. ■ 

Remark 10 Proposition fy) states that the probability density p^,. T is localized within the x 1 
interval r ■ [v\, V2] for r — ► ±00. In case of < k\ < k-2 it is localized in the half space x 1 < 
for t — > —oo and in the half space x 1 > for r — » 00. The localization is right moving. In case 
of k\ < k% < it is localized in the half space x 1 > for r — > —00 and in the half space x 1 < 
for t — > 00 . It is left moving. 
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4 Plane waves for V > 2k 

Let ip = exp(—iujx°)f(x 1 ) with / : R — > FT" be a single frequency solution of (||) with the 
frequency ui > k. With some a, /3 € C the equation ip = allk + PU-k holds on M_ := 
{p e M : x 1 (p) < 0}. Here fc > is determined by w(fc) = uj. On Af+ := {p e Af : x 1 ^) > 0} 
the function exp(iVx°)ip equals a single frequency solution of equation (gj) with A = 0. Its 
frequency reads a/ := uj — V. 

Klein's phenomenon occurs for uj' < —k. Thus we restrict our discussion to the case k < 
uj < V — k. This implies V > 2k. In that case 

ip = exp(-iVx°) tyV q + 5V- q ] (on M+) 

with 7, S G C and with q > being determined by —uj{g) = u>' = ZU(fc) — V. The constants 
a, (3, 7, 8 are restricted by the condition that ip is continuous, which is equivalent to 

aU k (0, 0) + pU-k(P, 0) = 7^(0, 0) + SV. q (0, 0). (3) 

Since Uk{0, 0) and C/_fc(0,0) are linearly independent, this system of linear equations for 
(a, /3, 7, <5) is of rank 2. Thus the space of single frequency solutions j/j to (||) is 2-dimensional. 
Within this space there are several Id subspaces of particular physical importance. 

One of these spaces comprises the single frequency solutions with 6 = 0. Its relevance 
emerges from the asymptotic behaviour in time of the wave packets formed from these solu- 
tions. As will be discussed in the next section, the x° = r restrictions of such wave packets 
localise for r — > — oo within the half line x 1 < 0. Thus these packets have a right moving 
incoming asymptotics. Similarly, the a — plane waves are the building blocks, from which 
incoming left moving packets are formed through superposition. Wave packets built from either 
a = or S = plane waves correspond to the "incoming" scattering solutions of the general 
quantum scattering theory. |l3| They have a well defined half space localization and direction 
of movement for r — > — oo. "Outgoing" scattering solutions approaching a wave packet which, 
for t — > oo, moves out exclusively towards x 1 — > — oo, are obtained from 7 = 0. These are the 
solutions, on which refs. H, || base their discussion of Klein's phenomenon. Finally, "Outgo- 
ing" solutions, which approach a wave packet, that, for r — > oo, exclusively moves out towards 
x 1 — > oo, are obtained from (3 = 0. 

Up to a constant factor the single frequency solutions with (5 = are given by the following 
lemma, which follows from the continuity condition (0) with a — 1 and (5 = 0. 

Lemma 11 Let V > 2k > and k > be such that ZU(fc) < V — k holds. Let the function 
ip : M — » W be continuous and such that for some q > and some r,(eC 



Uk + rU~k on M_ 

cxpt-iUa; )^ o?i Af+ 



holds. Then is a solution of (j^J i/ and onZj/ if ij, Sj, and 3j ZioZd. 
1) q is the unique solution ofu~j{k) +U(q) = V in M>o, 

3)t = t(k) : =-2^ffiM^. 
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Notation 12 The solution ip of which is given by lemma is denoted by U™ in what 
follows. 

Remark 13 Observe that due to V 2 — (k — q) 2 = 2(k 2 + ZJ(k)uJ(q) + kq) > and due to 
V + k — q = uj(k) + k + uj(q) — q > uj{q) ~ q > the inequalities r(k) < and t(k) < hold. 

Remark 14 The wave number q is given explicitly through q = s(k), with the differ entiable 
function 

s : (0, \jv 2 - 2nV) -> (0, \jV 2 - 2kV), k h-> ^Jv 2 - 2Vu:(k) + k 2 . 



The function s is a monotonically decreasing bijection. It has the fixed point = fy 1~ (y) ■ 
Thus if k < ko then q > ko and if k > ko then q < kg. These inequalities will show up in the 
group velocity of narrow wave packets through either an acceleration or a deceleration upon 
transition through the potential's singularity at x 1 = 0. 



Proposition 15 For the current vector field j V in 

TTj ur = t(k) 2 [ZJ{q)d Q + qdi] + 2nr(k)Q(-x 1 ) [cos^x 1 ) - l] d 

holds, jjjin is differentiable and div^^in) = everywhere on M. 

Proof. From remark (|^) one obtains by inserting the restrictions of f7£™ to M± into the 
current's definition the equations 

J [uJ(k) (l + r(k) 2 ) + 2nr(k) cos(fca; 1 )] do + k (l - r(k) 2 ) di on M- 
^ Jc/ r " { t(k) 2 {Lo{q)do + qd x } on M+ ' 

The continuity of U™ implies the continuity of j\J in , which in turn is equivalent to the equations 

U(k) (1 + r(fc) 2 ) + 2nr{k) = uj(q)t(k) 2 , 
k(l-r(k) 2 ) = qt(k) 2 . 

(They are easily checked by inserting the explicit expressions for r(k) and t(k).) From this the 
formula for jjjin follows on M. Though is not differentiable where x 1 = 0, its current field is 

diffcrentiable in every p £ M because of [cos(fca; 1 ) — l] (p) = and (<9 M [cos(kx 1 ) — l]) (p) = 
for x 1 (p) — 0. Obviously, div(j[/in) = holds on M + U M_. Since j V in is differentiable on M, 
it follows that div^'y*™ ) = on M. ■ 

Remark 16 From the j 1 -continuity condition (l — r(fc) 2 ) = qt ^ > and from r{k) < the 
bounds — 1 < r(k) < follow. 



Remark 17 Since the Lie bracket 
invariant. 



9o,ju irl = 0? the current vector field jy™ is x -shift 
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Definition 18 Let j be a differentiable vector field on M. Let j p : L — > M obey the differential 
equation j p = j o j p and the initial condition 7(0) = p. The open interval I Cl is assumed to 
be maximal. Since j is differentiable, -f p is unique. 7 p is called the maximal integral curve of j 
through p G M and the set J P (I) C M is called the orbit ofj p . If j is the current vector field 
of a solution ip of the Dirac equation, then the orbit of j p is called the Bohmian trajectory of 
if through p. 

Remark 19 jijin is a bounded and everywhere future directed time like vector field with pos- 
itive, constant component j^ in . Thus along a Bohmian trajectory of Ujp the function x 1 in- 
creases with x° , i.e. the trajectories move towards the potential step from x 1 < 0. To be more 
specific: On M + the equation j V i n — t(k) 2 ju q holds. On M_ the current j V in is the sum of 
t{k) 2 ju q and the non constant do-directed, bounded vector field 2 ^ L r(k) [cos(fca; 1 ) — l] do- Due 
to r(k) < the additional term belongs to R>o ■ do, which implies 

Ak 

< t(kff Uq < f ur < t(kff Uq + — \r(k)\ on M_. 
From this and j V vn = t(k) 2 ju q (on M + ) it follows that, due to 

GUuinJufn) > t(k) 2 G(j Uq ,J Uq ) = t(k) 2 ^ > 0, 

j/jir, is globally time like and future directed. Thus the current has no zeros and is nowhere 
light like. 

fan 

Remark 20 The velocity vector field v V i n := -^r—d\ of j V i n relative to the inertial frame 
(do, d\) is given by 

-1 

ft. 

On M + it is the constant field =^jdi and on M_ it oscillates between the positive bounds 

01 < Vrjin < ——01. 

From the bounds of u^.n it is obvious that a Bohmian trajectory of U™ cannot have a higher 
velocity within the range M_ than it has within M + , where its velocity is less than 1. An explicit 
formula for the Bohmian trajectories of U™ is given by the following proposition. 

Proposition 21 The Bohmian trajectory I\ T of U™ through po ■= (r, 0) is the set of all 

p G M, on which 

x o_ t= ^^i _ Mg^ij l kx i _ sin(jb .i)l (4) 
q ktyky 

holds. {I\ T } rgR is a disjoint covering of M and I\ T n {p € M : x 1 (p) = 0} = po- 



v ur = q 
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Proof. The current j := j V in — ado + bdi +c9(— — cos(fca; 1 ))9o with na := uj(q)t(k) 2 , 
nb := qt(k) 2 , ixc := —2nr(k) is differentiable on M. With 7 M := x^ 07 the differential equation 
7 = j o 7 decomposes into 7 1 = 6 and 

_ / a for 7 1 > 



7 I a + c(l - cos(fc7 1 ) for 7 1 < 

The initial condition 7 1 (0) = thus implies 7 1 (A) = bX for all A £ 1. Inserting this into the 
equation for the component 7 one obtains from 7 (0) = r 

,.% _ J aX + r for A > 

7 > ~ \ (a + c) A - ^ sinO&A) + r for A < 

Thus the maximal integral curve 7 Po of j through po reads 

7po : R - M, 7P o (A) = (r, 0) + A • (a, b) + 9(-A) (c • (A - !EjM), 
Its orbit is the set 1^ T of points p E M, on which 



x° - t = yx 1 + ^jOi-x 1 ) [kx 1 - sinikx 1 )} 
h kb 



holds. Obviously T^ T — (r,0) + Tj, holds. As x°\ r is expressed in terms of a; 1 L , the 
trajectories T}- iT and T^o therefore do not intersect for r 7^ 0. In particular, T^o intersects 
{p e M : x 1 (p) = 0} only at the single point (0,0). Because x 1 ] — > R is a bijection, the 
Bohmian trajectory of j through an arbitrary point p is obtained by the proper choice of r. 
Thus {I\ r } TgR is a disjoint covering, a fibration of M. m 

Figure (1) shows the trajectory Tk.o for V = 2.25k and k = § within the space time region 
where -600 < kx° < 200 and -50 < kx 1 < 50. 



5 Incoming localized solutions for V > 2k 

Notation 22 Let a : I — > C 6e continuous on a closed interval I = [fci, fa] C R>o o/ positive 
real numbers such that on I the inequality u7 < V — k holds. The constants V and k obey 
V > 2k. Then U in [a] :M->W denotes the wave packet U tn [a] := jj dp(k)a(k)U l k n . 

Proposition 23 U m [a] is continuous on M and differentiable on U . It is a solution of equa- 
tion (jl). Let a r : [—£2, — fa] — > C and at ■ [s(fa), s(fa)] - * C be defined by a r (k) :— r(—k)a(—k), 
a t (q) := %t(h)a(k) with k = ^/V 2 - 2VoJ(q) + q 2 . Then the equations U in [a] = U [a] + U [a r ] 
on M_ and U m [a] = cx_p(-iVx°)V [a t ] on M+ hold. 

Proof. From [/| n = U k + r{k)U- k (on M_) there follows on M_ 



' dp(k)a(k)Ur = ( 2 dp{k)a{k)U k + [ ' dfi(k)a(-k)r(-k)U k . 
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Figure 1: i\ for k — k/2 and V = 2.25k 



This proves U m [a] =U[a] + U [ a r ] on M_. 

Similarly, with q = s(k) := ^V 2 - 2VZJ(k) + 



k 2 there follows on M4 

„0\ 



fci 



dn(k)a(k)Ul n = exp(-iya; u ) / d^(k)a(k)t(k)V s{k) . 

Jki 



The function s is defined implicitly by 



V. This implies first s = s 1 and second 



cEZJ = — d (ui o s). Substitution of the integration variable k by q — s(k) then yields 



A-2 



dk 



fci 



2 ^a(fc)*(fc)V a(fc) 



s(fc 2 



du(q) q 



2q s(q) 
Q 



a(s(q))t(s(q))V q 



s(fcl) 
s(fci) 

d^{<l)^-Ta(s(q))t(s(q))V q . 

s(k 2 ) s Wj 



This proves C/ m [a] = exp(— iT^x )^ [a 4 ] on M + . The statements about the continuity and 
differentiability follow from the continuity of a r and a t through application of elementary 
analysis theorems. ■ 



Remark 24 The localization ofU m [a] for t — > ±00 is immediate from proposition fj23j). Since 
a r has its domain within R<o and a and at have their domains within M>o, the wave packet 
U [a r ] is left moving, while U [a] and U [at] are right moving. At large negative times U m [a] 
approximates U [a], which moves in through M_ towards x l — 0. At large positive times U m [a] 
approximates U [at] + U [a r ] . The transmitted wave U [at] moves away from x 1 = through M + 
and the reflected one, U [a r ], moves away from x 1 = through M_. 
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Rigourous analysis of whether the Bohmian trajectories of U m [a]| M may be connected 
with the trajectories of U m [o]\ M across x 1 = in a unique way is left to further investi- 
gation. As an indication that this should be possible, we present some trajectories, computed 
numerically. 

For V — 4k we choose for several values of the constants A > and K > the Fourier 
amplitude 

a : [k%, k 2 ] — > K, a(k) := exp 

with ki = K — 2A, k 2 = K + 2A. The domain of the function s of remark is < k < y^K 
and its fixed point is ko = \/3k. Clearly, the constants K and A have to be chosen such that 
[k\, k 2 ] is contained in the domain of s. First the wave packet U m [a] is computed numerically 
in that space-time region, where it hits the potential step. Second from the associated current 
field some Bohmian trajectories are computed by numerical integration. The starting points 
are chosen within S T around the center of localization of pu in [a],T- The initial time r is such 
that the main bulk of the probability distribution pu™\ a ],T nas n °t ve t arrived at the potential 
step. 

Figure (2) shows some trajectories for A = 0.1k and K — 0.3k within the space-time region 
where —180 < kx° < 200 and —100 < nx 1 < 300. Since the wave numbers k from the domain 
of a obey k < fco the domain of a is mapped into k > fco by s such that the transmitted packet 
is faster than the incoming one. This might be considered as what is left over from Klein's 
paradox. 





Figure 2: trajectories for A = 0.1k and K = 0.3k 

Figure (3) shows some trajectories for A = 0.05k and K = 2.7k within the space-time 
region where —100 < nx° < 160 and —300 < nx 1 < 200. Since the wave numbers k from the 



-{k-Kf 
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domain of a obey k > fco the transmitted packet is slower than the incoming one. Many more 
pictures for related situations are contained in ref. jig]. 



KX 



100 



-100 




100 k x 



Figure 3: trajectories for A = 0.05k and K — 2.7k 
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